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Abstract. We discuss neutral B s d meson mixing and 6-hadron lifetimes from the 
perspective of recent lattice calculations. In particular, we consider matrix elements 
which can be combined with measured A¥ Sj j to constrain \Vtd\, the lifetime ratios 
T(A b )/T(B%) and t(B~)/t(B%) and the lifetime difference, AT b JTb s , in the neutral 
B s meson system. 



| Based on talks given by the authors in the UK Phenomenology Workshop on Heavy Flavours and 
CP Violation, 17th - 22nd September 2000, Durham, England. 
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1. Introduction 



In this workshop, the lifetimes and mixings working group identified a 'wish 
list' of hadronic quantities whose accurate theoretical determination is crucial for 
understanding current and future 6-physics experimental results. We discuss several 
of these quantities from the perspective of recent lattice calculations: 

• matrix elements which can be combined with measured AM St d to constrain \Vtd\ 

• the lifetime ratios r(A 6 )/r(S2) and T ( B ~)/ T ( B d) 

• the lifetime difference, AT Bs /T Bs , in the neutral B s meson system. 

For a full survey of recent 6-physics results from the lattice see the reviews in [0-0. 



2. AM d and AM s /AM d 



The mass difference of the Bd-Bd system, AMd, constrains the poorly known CKM 
matrix element \V t d\- AMd has been experimentally measured to good accuracy. On 
the theory side, the main uncertainty comes from the long-distance strong-interaction 
effects in the matrix element 

M B M = (B d \QM\B d ), (1) 

which appears in the Standard Model prediction for AMd to leading order in an 
expansion in 1/M W [0J§. In Eq. (fj), Q Ld is the four-quark operator [67^(1 — 7 5 )<i] [^7^(1 — 
7 5 )<i] and /i is the renormalisation scale. 

An alternative approach, in which many theoretical uncertainties cancel, is to 
consider the ratio, AM S / AMd, where AM S is the mass difference in the neutral B s — B s 
system. In the Standard Model, one has 
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where Ql 3 is the same operator as Qi d with d replaced by s and where the 
renormalisation-scale dependence of these operators cancels in the ratio. Because the 
unitarity of the CKM matrix implies |Vt s |^|Vd,| and because \V c b\ can be accurately 
obtained from semileptonic B to charm decays, a measurement of AM S / AMd determines 



td\ 



This is experimentally very challenging because of the rapid oscillations in the 
B Q S —B° S system. Nevertheless, the experimental lower bounds on AM s /AMd already 
yield interesting constraints on the b—d unitarity triangle 0. 

The matrix elements in Eqs. (|l|) and (Q) are traditionally parameterised by 



M B U) = {B q \Q L U)\B q ) = -M&IBbM 



(3) 



with q = d or s, where the i?-parameter, B Bq , measures deviations from vacuum 
saturation, corresponding to B Bq — 1, and f Bq is the leptonic decay constant. One also 
usually introduces a renormalisation-group invariant and scheme-independent parameter 
B Bq , which to NLO in QCD is given by 



Bt° = CMB B U), 
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where Cb is the two-loop Wilson coefficient calculated in the same scheme as the matrix 
element f7|,|T0|. 

Because Mb 3 and Ms d are measured experimentally, one needs to calculate £ and 



fB d \jB B x ° non-perturbatively to determine \V t d\ from experimental results for AM S / AMj 
and AM d . 

There are three recent quenched lattice calculations of £ |[ The APE |13|] and 
UKQCD (with preliminary results reported previously in fl4}| ) Collaborations use 
relativistic formulations of quarks and obtain £ at the physical B meson mass by 
extrapolating from heavy- meson masses around that of the D. Gimenez and Martinelli 
(denoted as GM in the following) [15| calculate this quantity in the static limit where 
the b quark mass is taken to infinity. The results from these three groups are: 

f 1.16(7) (APE) 

1.16(2)11 (UKQCD) (5) 
[ 1.17(3) (GM), 

where the first error is statistical and the second is systematic. In order to avoid the 
situation where discretisation errors, which can be significant in lattice calculations 
involving propagating heavy quarks, are out of control, both APE and UKQCD perform 
numerical simulations with several meson masses straddling the D and then extrapolate 
to the B meson mass using HQET-inspired relations. APE use a fully 0(a) (a is the 
lattice spacing) improved fermion action |1| in which the leading discretisation error 
is 0(a 2 ). UKQCD use a mean-field improved fermion action Jl7|], and the leading 
discretisation error here is formally 0(aa s ), although it might be numerically smaller. 
However, it should be noted that both groups have 0(aa s ) discretisation errors in £, 
because the four-quark operators Qi q are not fully 0(a) improved in these calculations. 
These 0(aa s ) errors are absent in the decay constants, and hence fB s /fB d , as calculated 
by APE, because they use the improved lattice axial current in the calculation. 

In Eq. (^|), UKQCD's result for £ has a much smaller statistical error than that 
obtained by APE, although these two groups have very similar statistics in the Monte- 
Carlo simulations. This is because UKQCD calculate £ by performing heavy-quark- 



mass extrapolations in the ratios fs s y/ M Bs / fB d J M Bd and B Bs /B Bd , in which the 1/ttiq 



corrections cancel significantly, especially for the decay constants f!2] , |14 |. This procedure 



determines £ to a better statistical accuracy than that calculated by extrapolating 



fB s V M B s , fB d \/M Bd , B Bs , and B Bd individually. 

The above three studies suggest that £ has small systematic uncertainties arising 
from discretisation effects and heavy-quark-mass extrapolation^], within the quenched 
approximation. The UKQCD result includes a systematic errorffl not estimated by APE 



§ The SU{3) breaking ratio M.b b / M-B d has also been studied in ]TlL[12[ by using a different method, 
in which one does not need to calculate £ . However, results from this method have large uncertainties. 
|| As mentioned above, the authors of (l2) look at the 1/toq corrections for (fs s \/ Mb, ) / \J Mb a ) 
and BB s /BB d and find that they are very small for both quantities. 

% This systematic error looks small on £. However, it should be noted that it is the error in £'s deviation 
from unity. 
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and GM. This is due to the uncertainty in the determination of the lattice spacing (see 
below for details), which introduces a variation of the strange-quark mass. 

Quenched Chiral Perturbation Theory (q%PT) predicts that quenching errors in 
BbJ ' Bs d are small if the couplings in the theory are constrained by large- N c arguments 
fT"8|| . Recent numerical studies with two flavours of dynamical quarks show little 
variation in fB s /fB d compared to its quenched value |T^|-f2^]. Therefore, quenching 
effects should not be significant in £. For more detailed discussions on this issue, please 



refer to [23 



APE |13[] and UKQCD [0 also calculate fB d yB B ° in quenched approximation. 
The results are: 

206(28) (7) MeV (APE) 



fB d ^B B lo d | 211(21)(28) MeV (UKQCD), ^ 

where the first error is statistical and the second is syste matic. 

The systematic error in APE's result for fg d J B B ° reflects the typical size of the 
0(a 2 ) discretisation effects in fg , a very small (~ 2%) uncertainty in the inverse lattice 
spacing and the uncertainty in matching the lattice-regularised four-quark operators 
onto the NDR-MS scheme, with these three errors taken in quadrature. 

In addition to the 0(a 2 ) and 0(aa s ) discretisation errors in f Bd , UKQCD 
investigate the uncertainties in the procedure of operator matching, heavy-quark-mass 
extrapolations and a ±7% uncertainty in the inverse lattice spacing. The systematic 
error in UKQCD's result is obtained by taking these uncertainties in quadrature. 



The overall systematic error on fB d yB B l d obtained by UKQCD is much larger 
than that obtained by APE (see Eq. @), mainly because of the estimate of the error 
in the determination of the inverse lattice spacing, a -1 . APE use a -1 set by the 
method described in |24j] for their central value, and the ones set by M p and to 
get the systematic uncertainty. This method results in a 2% systematic error. UKQCD 
determine a -1 in conjunction with the strange-quark mass by fx and Mk, and then 
vary a" 1 by ±7%, a range which covers the typical variations of a -1 set by gluonic 



and light-hadron spectral quantities with the same action used in their calculation [25 



This procedure introduces a 9% effect on their result for fB d yB B °. In addition, the 
systematic uncertainty from the heavy-quark-mass extrapolations (~ 3%) which is not 
estimated by APE, and the 0(aa s ) discretisation errors in fs d (~ 5%) which are not 
present in the APE calculation, further enlarge the gap. 

Both APE and UKQCD do not attempt to quantify the discretisation errors in 
the B parameters, as this requires one to consider the mixing with dimension-seven 
four-fermion operators and is therefore very involved. Nevertheless, B Bq are ratios 
between closely related matrix elements, hence one expects that these errors might 
cancel significantly. As mentioned above, the fermion actions used by APE and UKQCD 
have different discretisation errors. However, as displayed in Table |I| and Figure [I], 
their results of the B parameters show very good agreement. Furthermore, UKQCD 
perform the calculations at two lattice spacings and observe that the change in Bb c 
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Table 1. Comparison of the latest lattice results for Bb 9 parameters renormalised in 
NDR-MS scheme at m&. The first error on each result is statistical, while the second 
one is systematic. For a more complete list of these results, including the earlier 
calculations, please refer to 

B Bd (m b ) B Bs (m b ) 



APE (a" 1 
UKQCD [a 
UKQCD la- 
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NRQCD (a^ 1 
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Figure 1. Comparison between APE and UKQCD 's results of renormalisation group 
invariant B parameters at different meson masses (from p"3|). The inverse lattice 
spacing in this plot is approximately 2.7 GeV. Mp is the heavy-meson mass in lattice 
units (i.e., one should multiply it by 2.7 to obtain the mass in GeV in this case). The 
open symbols are the results from numerical simulations, and the closed symbols are 
those extrapolated to the Bd meson mass. 



is insignificant. Finally, both groups obtain B Bq compatible with those in a slightly 
earlier work by Bernard, Blum and Soni (BBS) who use a less improved fermion 
action but extrapolate their results to the continuum limit (See Table 0). Recently, 
these B parameters have also been calculated using NRQCD to describe the heavy 
quarks [2£|p7|, with the heavy-quark masses in numerical simulations around that of 
the b. The central values of these results are 2 standard-deviations lower than those 
obtained by APE and UKQCD. However, as Table shows, when the systematic errors 
are taken into account, all these results are compatible. 
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UKQCD also normalise fB d yB B l ° by fo a because some systematic errors cancel, 
and there have been preliminary experimental results of fo s P§|, 2~9|j. They obtain [12] : 



^^ = 0.89(7)^, (7) 

JD S 

where the first error is statistical and the second one is systematic. 

While quenching errors in Bb 9 are predicted to be within a few percent by qxPT 
with reasonable ranges of the couplings in the theory [Tf8| , they can be significant for the 
decay constants. Recent numerical studies with two flavours of dynamical quarks flT9|-f22] 



show visible increases on the decay constants. Again, please refer to for more 
detailed discussions on quenching errors. 

Finally, there are on-going analyses J]],[30],|3l| in which results from the static-limit 



and relativistic-heavy-quark calculations are combined. In such analyses, fs d y B B ° and 
are obtained by interpolations between the charm- mass region and the limit of 
infinite heavy-quark mass. 

3. Spectator Effects on 6-Hadron Lifetimes 

The combination of operator product and heavy quark expansions predicts that the 
ratio of lifetimes of two hadrons, Hi and H 2 , each containing a single 6-quark, is given 



by 32 



r(ffo _ i , Am - Am vim - Am , 

Here the leading 1 on the right hand side arises from a universal first term describing 
free 6-quark decay. There is no term of order l/m& in the expansion. At order l/ml, the 
A and A terms are matrix elements in the heavy quark effective theory of the kinetic 
energy and chromomagnetic moment operators respectively between the corresponding 
hadron states. They can be fixed from hadron mass formulas and mass splittings. The 
coefficient cq is calculated in p3|-|3"6|j . Therefore, for the cases of the ratios of B~ and B® 



lifetimes or A b and B^ lifetimes, we have calculated values to compare to experimental 
results: 

t(B-) (l+0(l/m 3 b ) r(A 6 ) ( 0.98 + 0(1 /m 3 b 



r(B° d ) 1 1.066(20) [37] r(B° d ) 1 0.794(53) 



(9) 



For the 5-meson ratio, there is no difficulty However, for the A^ to B ratio, the 0(l/mf) 
terms apparently must account for almost all of the deviation from one. 

Some work has addressed the question of whether "spectator effects" in the 0(l/ml) 
terms can be responsible for such large deviations. In the operator product expansion, 
the effects of the diagrams in figure |2| lead to the appearance of AB = four-quark 
operators of the form bTiqqVib. These diagrams are the first time the effects of spectator 
quarks enter explicitly. Moreover, the diagrams are one-loop, compared to the two-loop 




Figure 2. Spectator contributions: diagrams like those on the left, where the box 
denotes a |A£>| = 1 transition, produce 4-quark operators in the operator product 
expansion. I\ denotes some combination of Dirac and colour matrices. See 132]. 



diagrams which generate 0(l/m%) terms in the expansion. There is thus a hope that 
these spectator effects could be large. The four quark operators are classified as 

0\ = b Llfl q L q L Yb L Tl = b Llil T a q L q L YT a b L (10) 
0\ = b R q L q L b R T| = b R T a q L q L T a b R (11) 

where q = u,d, s. The spectator contribution to the decay rate of a hadron H depends 
on (H\O spec \H), where 

O spcc = F ul (z)0^ + G ul (z)T? + Yl ( F i<,( z )°i + Gn{z)T?) . (12) 

i=l,2;q=d,s 

Here z = m 2 c jm\ and the F's and G"s are coefficient functions known at leading 



order [32 



Given the coefficient functions, the remaining ingredient is the evaluation of the 
four quark matrix elements. Lattice QCD simulations of these [^8| -|40|| are reported here 
(for sum rule calculations see f4T,[4"2"f). To date, these calculations are still exploratory 
for the A b meson, but the results are encouraging and the calculations deserve repeating. 

It is convenient to parameterise the matrix elements. For mesons we use factors Bi 
and e i7 

1 -(B\Of\B) = l^B t , J-(B\Tf\B) = S^e t , (13) 



2m R 1 8 2m R o 

for i = 1,2. From vacuum saturation or large iV c arguments the expectation is that 
Bi ~ 1 and the are small. For the A& baryon, heavy quark symmetry implies that 
only two matrix elements need to be considered (up to l/m& corrections) ||32|| . The 
corresponding parameters, L\ and L2 are defined byQ, 

1 (A 6 |0?|A 6 ) = ^L 1; -^-(A 6 |Tf |A 6 ) = ^L 2 . (14) 



2m b 8 2m b 

The conversion to the parameters used in J|^] is B = — 6L1 and r = —2L 2 /Li — 1/3 
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destroy B 0?,lf created destroy A b Of,T? create A b 

Figure 3. Lattice quark propagator contractions needed to compute spectator effects 
on B and lifetimes. 



In terms of these parameters the lifetime ratios can be expressed as, 
r(B-) 



r(B d 



^ — 1 + a x ex + a 2 e 2 + a 3 B 3 + a A B 4 , (15) 
■ = 0.98 + fexei + b 2 e 2 + 6 3 ^i + hL 4 , (16) 



r(B° d ) 

where tiny B 12 terms are neglected in the second equation. The coefficients a { and bi 
are known perturbatively, 

,697 6i = -0.175 

1 QF> h„ = n 1 Q£ 

(17) 



These values are quoted in the MS scheme at a scale // = mj and are to be combined 
with matrix elements evaluated in the same scheme. 



CLi = 


-0.697 


h = 


-0.175 


a 2 = 


0.195 


b 2 = 


0.195 


a 3 = 


0.020 


h = 


0.030 


04 = 


0.004 


b 4 = 


-0.252 



For the B meson, the lattice calculations |38| have been done in the quenched 
approximation with a static b quark at an inverse lattice spacing of a -1 = 2.9 GeV 
(corresponding to an input lattice coupling parameter (3 = 6.2). The results have been 
extrapolated from the light quark masses actually simulated (for technical reasons, one 
cannot simulate with realistically light quark masses) to the chiral limit. For the 
meson, the matrix element involves a more complicated set of lattice quark propagator 
contractions, as shown in figure |3|, and the calculation is still exploratory |39|, flOf . It is 
also done with a static b quark, but on a coarser lattice (a -1 = 1.1 GeV or (3 = 5.7) 
and using a stochastic technique to calculate the light quark propagators |43|] . No chiral 
limit is taken in this case, the results are quoted for pion masses given by am w = 0.52(3) 
(case A) and am n = 0.74(4) (case B). The results are: 

B x = 1.06(8) ei = -0.01(3) 

B 2 = 1.01(6) e 2 = -0.02(2) 



(18) 



for the B and 



, -0.31(3) A f 0.23(2) A 

Li = { 1 ^ L 2 = I 1 ; (19) 

1 -0.22(4) B I 0.17(2) B 
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for the A;,. Combining these with the calculated coefficients and bi leads to, 

t(B~) T(Ah) [ 0.91(1) A 

^J- = L03(2)(3), = 1 K } (20) 

r(B° d ) r(B° d ) \ 0.93(1) B 

We see that the answer is close to unity as expected for the B meson ratio. For the 
Afe to B ratio, about 40% of the deviation from unity required to match experiment 
is reproduced. Given the exploratory nature of the A& calculation, this indicates that 
l/mf spectator effects could be large enough to suppress the A& lifetime compared to the 
B meson lifetime and further lattice studies are warranted. Moreover, next-to-leading 
order calculations of the coefficient functions are in progress [ 44j . Experience with the 
B s lifetime difference, where the next-to- leading effects are significant, provides further 
motivation. 



4. Width Difference of B s Mesons 

The decay width difference of B s and B s mesons arises from the forward off-diagonal 
matrix element of a time ordered product of weak effective Hamiltonians mediating 
b quark decay. Applying the heavy quark expansion leads to an expression of the 
form 



2 

Ar ^ = ioC 6 I V <* V ™\ 2 (G(z)(QL s (m b )) - G s (z)(Q s (m b )) + 5 l/m Vl^A~z) (21) 
12ttM Bs v ' 

where z = rri^/ml and angle brackets denote matrix elements between B s and B s , 
(O) = (B s \0\B s ). The four-quark AB = 2 operators Ql s and Qs arise at leading order 
in the heavy quark expansion, while 8u m denotes l/m^ corrections involving further 
operator matrix elements. The coefficients G and Gs are known at NLO in QCD ||46|| . 



Ql s is the same single operator that contributes to the mass difference, AMe a , while 
Qs = (bs)s- P (bs)s-p. (22) 

Notice that (Ql s ) is actually the matrix element M.b s defined in Eq. @. Using the 
standard parameterisation, 

k y\f2 

(Qs(jt)) = -~flM 2 B -2* Bsfji), (23) 
3 (mj + m s y 

and Eq. (|), one can write, 

Ar Bs G\m\ 2 2 

\VcbVcs\ T Bs f Bg M Ba 



B. 



12tt 



2 \ 

^ G{z)B Bs + u - G s {z)—^—B s + SymVT^Iz' . (24) 



5 r ( \ M b 

„ y „ Ds , -G S {Z)- 1 ; xn -PS-l-Ol/ 

\3 3 [Tni ) + m s y / 

The coefficients Gig\ and matrix element parameters B Bs and Bs depend on a 
renormalisation scale \i such that the result is in principle \i independent. Results 
quoted below have \i set to the 6-quark pole-mass, mj, taking m& = 4.6 GeV. Note that 
different authors make different choices for the definitions of the quark masses appearing 
in the defining equation for Bs- Here they are the pole masses. 
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Input from lattice QCD is needed for the values of f Bs and the parameters B Bs and 
B$- From the lattice viewpoint, it is convenient to trade uncertainty in determining 
fs t |||23| for the appearance of extra CKM factors by considering Ar_B s /AMe s . The 
mass difference is proportional to (Ql s ), so that the ratio AT B J AM Ba depends on the 
quantity 

(Qs) 5 B s (m b ) 



TZ(m b ) 



Ml 



(25) 



(Ql s ) 8 B Bs (m b ) (m b + m s ) 2 ' 
Systematic errors from uncertainty in determining the lattice spacing and from 
quenching effects should be reduced in this dimensionless ratio of similar matrix 
elements. In the absence of an experimental measurement for AM Bs , one further uses 
the ratio £, as defined in Eq. (|2|), which is quite well determined by lattice calculations 
(see section 0). In this way one arrives at an expression [|? 

Ar Bs 4vr M Bs m\ V cb V cs 2 r Bs AM Bd 2 , ~ 

— — — — C {G(z)-G s (z)n(m b ) + 5 1/mi 



r. 



3 M Bi m 2 w 



V tb V td 



(26) 



r) B (m b )S (x t ) 

The quantities rj(m b ) and So(xt), where x t = mf/m^, come from the expression for 
AM Bs and are known factors 0,0. At leading order in l/m b the non-perturbative 
contribution is isolated in 1Z(m b ). At order l/m b , one needs both the explicit 8u m piece 
together with the implicit m b dependence from the matrix elements of Ql s ,s in the ratio 

n. 

Three groups have recent calculations giving values for 1Z. They use different lattice 
formalisms, but give very consistent results for this dimensionless ratio. Becirevic et 
al. [|7 use heavy quarks at around the charm mass and extrapolate to the h. The 
Hiroshima-KEK group |26[] use lattice NRQCD, simulating directly at the b quark mass. 
Finally, Gimenez and Reyes [pC, 49] work in the lattice static quark theory, and hence 
include 1 / m b corrections to the matrix elements as a systematic error. All groups have 
results from quenched simulations, but Gimenez and Reyes 0] also have results with 
two flavours of degenerate sea quarks. 



H(m b ) 



[ -0-93(3)(?) 
-0.91(5)(17) 
-0.95(7)(9) 
-0.97(5)(15) 



extrap c —* b [47] 
NRQCD |2§ 
static 6 

static b, rif = 2 



(27) 



The variation of about ±3% in the central values for TZ is smaller than the corresponding 
uncertainty of about ±10% in the values of B Bs and B$- 

Since the ratio G(z)/Gs{z) is 3-4% []4"6[1 , the contribution from (Qs) dominates 
that from (Ql s ) in A.T Bb /T Bs . However, the l/m b corrections, currently estimated 
using factorisation, cancel significantly against the (Qs) term []50[| , so that two terms of 
order 0.1-0.15 combine to give a prediction AT Bs /T Bs ~ 0.05 from equation Eq. (pG|) 
with a large uncertainty. Moreover, using results for B^b s ,s) and f Bs in Eq. (|24l) instead 
of Eq. (I 



]) gives AT Bs /T Bs closer to 0.1 [50]: here there is extra uncertainty from 
the lattice determination of f Bs , although the 5i/ m uncertainty dominates. To make 
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progress, better input for the matrix elements in S 1 / m is vital, together with better 
knowledge of fs s '■ both may be addressed by future lattice calculations. Note also that 
the 'ratio' form of Eq. ( p6|) assumes that AMb b (and AM^J are given by their standard 
model expressions. Ar Bs by itself is dominated by tree level physics and so is expected 
to be less sensitive to new physics than the mass differences. This would favour working 
with the expression in Eq. (^) once fs s is better known. 
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